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Abstract 

A general class of cosmological models driven by a nonlocal scalar field inspired by 
the string field theory is studied. Using the fact that the considering linear nonlocal 
model is equivalent to an infinite number of local models we have found an exact special 
solution of the nonlocal Friedmann equations. This solution describes a monotonically 
increasing Universe with the phantom dark energy. 

1 Introduction 

Recently string theory and brane cosmology have been intensively discussed as promising 
candidates for the theoretical explanation of the obtained experimental data (see for example 

mm)- 

The purpose of this paper is to present new results concerning studies of nonlocal linear 
models in the Friedmann-Robertson- Walker Universe. These models are inspired by the 
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string field theory (SFT) (for review of the SFT see [7]). A distinguished feature of nonlocal 
linear and nonlinear models [8] [21] is the presence of infinite number of higher derivative 
terms (note also nonlocal models in the Minkowski space-time [22] [29]). For special val- 
ues of the parameters these models describe linear approximations to the cubic bosonic or 
nonBPS fermionic SFT nonlocal tachyon models, p-adic string models or the models with 
the invariance of the action under the shift of the dilaton field to a constant. The NonBPS 
fermionic string field tachyon nonlocal model has been considered as a candidate for the dark 
energy [8]. 

Present cosmological observations [30] do not exclude an evolving dark energy (DE) state 
parameter w, whose current value can be less than —1, that means the violation of the null 
energy condition (NEC) (see [311 E2] for a review of the DE problems and [33J for a search 
for a super-acceleration phase of the Universe). 

Field theories, which violate the NEC [34, 35J, are of interest not only for the construc- 
tion of cosmological dark energy models with the state parameter w < — 1, but also for the 
solution of the cosmological singularity problem. A possible way to avoid cosmological singu- 
larities consists of dealing with nonsingular bouncing cosmological solutions. In this scenario 
the Universe contracts before the bounce \2\. Such models have strong coupling and higher- 
order string corrections are inevitable. It is important to construct nonsingular bouncing 
cosmological solutions in order to make a concrete prediction of bouncing cosmology. 

A simple possibility to violate the NEC is just to deal with a phantom field. In the present 
paper we consider nonlocal models which are linear and admit solutions, which are linear 
combinations of local fields. Some of these local fields are phantoms. Namely due to the 
presence of these ghost excitations such nonlocal models present an interest for cosmology. 

At the same time there are well known problems with instability of quantum models with 
phantoms, namely a lost of unitarity and so on. We believe that nonlocal SFT models in 
true vacua are stable with respect to quantum fluctuations. This question has to be consider 
in the full SFT framework and demands further investigations. We also believe that due 
to these string theory origin the corresponding nonlocal cosmological models, which are 
nonlinear in matter fields, have no problem with instability in the quantum case. In this 
paper we consider only the classical case and models, which are linear in a nonlocal scalar 
field. 

In our previous paper [15] as well as in paper [16] nonlocal linear models already have 
been studied. In [15] nonlocal linear model has been studied in the flat space-time and 
we have proposed special deformations of the potential, which allows us to get the same 
scalar field solutions in flat and nonflat (the FRW metric) cases. As result we have obtained 
nonlinear models in the FRW metric. In [T6] few exact solutions to linear model in the 
FRW metric have been found. In this paper we present a systematic method that permits 
us to transform the initial nonlocal system into infinity set of local systems. The choice of 
a local system is equivalent to the choice of a special solution of the nonlocal system. This 
approach allows us to use the standard method of analysis of the differential equations and 
in particular to find exact solutions. 

The paper is organized as follows. In Section 2 we describe string inspired models with 
quadratic nonlocal potentials. In Section 3 we assume that the metric is given and consider 
the equation of motion as an equation for the nonlocal scalar field. We construct solutions, 
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using eigenfunctions of the D^-operator with eigenvalues, belonging to the set of roots of the 
characteristic equation. In Section 4 we find values of the energy density and pressure for 
these solutions. In Section 5 we consider the Friedmann-Robertson- Walker Universe and find 
local models, which correspond to particular solutions of the initial nonlocal model. In the 
case of dilaton massless scalar field we construct the general solutions for the corresponding 
local model, which are the special exact solution for the initial nonlocal model as well. We 
analyze cosmological properties of the obtained solutions. 



2 Nonlocal linear models 

In this paper we consider a model of gravity coupling with a nonlocal scalar field, which 
induced by string field theory 

S = |^v^(^ + ^Q0F(-n 9 /M s 2 )0-A'^, (1) 

where g^ v is the metric tensor (we use the signature (— , +, +, +) ), D g = -^=d^—gg^ u d u , 
M p is a mass Planck, M s is a characteristic string scale related with the string tension 
a': M s = \j\fo! ', is a dimensionless scalar field, g^ is a dimensionless four dimensional 
effective coupling constant related with the ten dimensional string coupling constant go and 
the compactification scale. A = -^-A is an effective four dimensional cosmological constant. 

The form of the function F is inspired by a nonlocal action appeared in the string field 
theory. We consider the case 

F(z) = -£ 2 z+ 1 -ce~ 2z , (2) 

where £ is a real parameter and c is a positive constant. Using dimensionless space-time 
variables and a rescaling we can rewrite ([T]) for F given by (j2j) as follows 



S 

where 



d'x^g- M^R + |> D g( f> +\(<P 2 -c $ 2 ) - A') , (3) 



and nip = g^M^jM 2 s . Generally speaking the string scale does not coincide with the Plank 
mass. That gives a possibility to get a realistic value of A. 

The form of the term (e n9 0) 2 is analogous to the form of the interaction term for the 
tachyon field in the SFT action. The case of the open cubic superstring field theory tachyon 

corresponds to £ 2 = - 1/ (41n (^75)) ~ 0.9556 and c = 3 (see |25]-[2T|). 
The equation of motion for the scalar field has the following form 

(£ 2 D 9 + l)e- 2D ^ = c0. (4) 

The energy-momentum tensor 

T 2 6S ™ 
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has the following explicit form 



1,2 C 



-g a pc / dp [(e( 1+ ") n ^)(D s e( 1 -' , ) n ^) + (9 M e^ D ^)(^e( 1 -'') n ^)] + 



+ 2c J dp (d a e {1+p)U * cf) (d, 



Note that the energy-momentum tensor T a/ g includes the nonlocal terms, so the Einstein's 
equations are nonlocal ones. 



3 Generalization of Flat Dynamics 

3.1 Flat Dynamics 

In the flat case action ([T]) has the following form: 

Sfiat = \J d 4 x<pF(-n)<p. (6) 

If the scalar field depends only on time, then equation of motion (j3J) is reduced to the 
following linear equation: 

mm = o. (7) 

A plane wave = e at is a solution of (JTj) if a is a root of the characteristic equation 

F(a 2 ) = 0. (8) 
For a case of F given by equation (j7j) has the following form 

-£ 2 <9 o 2 + 0-ce- 2a o0 = O. (9) 

This equation has been analysed in detail in our paper [15] . Using the explicit form the 
function <p(t) we have found the solutions of equations of motion and the corresponding 
values of the energy density and pressure. In this paper we generalise these calculations for 
non-flat case. 

3.2 The equation of motion in an arbitrary metric 

Let us consider eq. flU). Really this equation is a consequence of the Einstein's equations, 
hence, both the metric and the scalar field <ft are unknown. We assume that the metric 
g^ u is given and consider eq. (jlj) as an equation in 0. 
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In this paper we study solutions in the following form: 

N 

= ^>n, (10) 

71=1 

where N is a natural number, n is a solution of the following equation: 

□ 5 0n = - «^0n, (11) 

and a n are solutions to the corresponding characteristic equation: 

F(a 2 n ) = -ea? n + l-ce-** = 0. (12) 

Without loss of generality we can assume that for any n and k ^ n the conditions o? n ^ a\ 
are satisfied. Indeed, if the sum ( TTUj) includes two summands 0^ and 4>k 2 , which correspond 
to one and the same a 2 ., then we can consider them as one summand = 0^ + 0fc 2 , which 
corresponds to a\. 

We start with construction of a solution to equation (@J in the case N — 1: 

□ 9 0= -a 2 0, (13) 

where a is a root of ( |T2l) . Note that this ansatz is widely used in studying of nonlocal linear 
models [6j [13l [16j [13, [181 [36]. Equation ( |T2l) has the following solutions 




«n= ±^/4 + 2£Wj , n = 0,±l,±2,... (14) 

where W n is the n-s branch of the Lambert function satisfying a relation W(z)e w ^ = z. 
The Lambert function is a multivalued function, so eq. f[T2l has an infinite number of roots. 
Parameters £ and c are real, therefore if a n is a root of f[T2l . then the adjoined number a* 
is a root as well. Note that if ot n is a root of (|T2l) . then -a n is a root too. 

If a 2 = «q is a multiple root, then at this point F(a^) = and -F'( a o) = 0- These 
equations give that 

2 1 1 

hence a 2 , is a real number and all multiple roots of F(a^) = are either real or pure 
imaginary. Double roots exist if and only if 

c= | e (2/^-l) (16) 

Note that the existence of double roots means that there exist solutions of equation (jlj), 
which does not satisfy of equation ([TBI , but satisfy the following equation 

(□ 9 + a 2 )(D 9 + a 2 )0 = O. (17) 



In the flat case an example of such solution is the function 0(t) = texp(at) (see 
All roots for any £ and c are no more than double degenerated, because F"(afy ^ 0. In this 
paper we consider such values of £ and c that equality (jTBl) is not satisfied and all roots are 
simple ones. Under this assumption we can consider the set of the solutions (jTOj) as a quite 
general solution. 
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3.3 Real Roots of the Characteristic Equation 



For some values of the parameters £ and c eq. ([12]) has real roots. To mark out real values 
of a we will denote real a as m: m = a. 

To determine values of the parameters at which eq. (fl2i) has real roots we rewrite this 
equation in the following form: 



£ 2 = g(m 2 , c), where g(m 2 ,c) 



e 2m 2 _ c 
m 2 e 2m 2 



The dependence of g(m, c) on m for different c is presented in Fig. (JTJ). This function has a 



maximum at rr? max 



provided c is such that W- 



m 2 = W i 



< —1, in the other words < c < 1. 



(19) 






Figure 1: The dependence of the function g(m 2 ,c), which is equal to ( 2 , on m at c = 1/2 
(left), c = 1 (center) and c = 2 (right). 

There are three different cases (see Fig. [T]). 

• If c < 1, then eq. f[T2"j) has two simple real roots: m = ±mi for any values £. 

• If c = 1, then eq. (IT2|) has a zero root. Nonzero real roots exist if and only if £ 2 < 2. 

• If c> 1, then eq. (JED has 

— no real roots for £ 2 > £ 2 M:E , where 

■ > 



Si/ 



1 — ce 2m ™ 



(20) 



— two real double roots m = ±m max for £ 2 = ^ ax 

— four real simple roots for £ 2 < ^ ax . In this case we have the following restriction 
on real roots: m 2 > | In c. 

Note that values of roots do not depend on H(t) and, therefore, coincide with roots in 
the flat case, which have been found in |15j . 
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4 Energy Density and Pressure 
4.1 General Formula 

Let us calculate the energy density and the pressure for the solution (fit)]) . Up to this moment 
we do not put any restrictions on the metric tensor g^, now we start to consider the case of 
the spatially flat Friedmann-Robert son- Walker Universe: 

ds 2 = - dt 2 + a 2 (t) (dxj + dx\ + dxf) (21) 

and spatially homogeneous solutions 4>(t). In this case 

T aP = g a(3 diag{£, -V, -V, -V}, (22) 

where the energy density S and pressure V are as follows 

S = 8k + S p + £ n l2 + £nll + A', V = Sk — S p + S n l2 — S n ii — A'. (23) 

Nonlocal term S n n plays a role of an extra potential term and £ n ii2 plays a role of an extra 
kinetic term. The explicit form of the terms in the R.H.S. of ( 123]) is [2H [29] as follows 

£k = |W) 2 , 



i 



Smi = c j (^+^0) ( - Ve^<P) dp, 
o 

i 

£ n i2 = - cj {de^4>) (de^<j>) dp, 



(24) 



where 



V = -d 2 - 3H(t)do, H =^. (25) 



For N = 1 we obtain 



£ = Efa) + A' = ^ ((doh) 2 ~ «M) + A', (26) 

V = P(0i) " A' = ^ ((<9 o 0i) 2 + afrt) - A', (27) 

where for arbitrary a 

r }a = Z 2 + 2?a?-2. (28) 

Note that considering the flat space-time [JS], we have introduced the parameter p a = 
a 2 rj a . The use of parameter r] a instead of p a is more convenient, because we do not need to 
consider the case a = separately. 
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Hereafter we denote the energy density and pressure of function <f>(t) as the functionals 
E(<p) and P(<fi), respectively. 

For the solution <p(t) = 4>i(t) + foit) it is convenient to write the energy density in the 
following form 

S = + <f> 2 ) +A' = + E{<f> 2 ) + E cross (<Pi, <h) + A', 

where the functional E cross (<pi, 2 ) is defined as follows: 

Ecross((pi, 02 ) = E kcr + E n i 2cr + E Pcr + E nilcr , (29) 

£ fccr = e 2 5o0i«9 o 2 , E Pcr = - 0!0 2 + ce--!-^^, 
i 



E nllcr =-cJ [(e^fc) V (e^<j> 2 ) + (e^fo) V (e^fa)} dp, 



o 

i 

E n l2 cr = — c 


Using ffl2l . we calculate E ni2 



dp. 



c e 



^ni2 cr = " 2 _ 2 ~d (f)id (f)2 = ~ £ <9 O 0l<9 O 02. (30) 

So 

^n/2„. + E kcr = 0. (31) 

The straightforward calculation also gives that 



Ka CT = - c e-^-^fafa + -± - 2 - 2 ^0i0 2 = - E Pcr . (32) 

«2 — Oi{ 

Therefore, we obtain that 

Ecrass (01, 2 ) = 0, and -P C r OSS (01,</>2) = 0, (33) 

where 

-P C ross(0i, 02) = E kcr + E n i 2cr — E Pcr — E n i lcr . (34) 

So, 

£(0i + 2 ) =E{<j> 1 ) + E{(j) 2 ), (35) 
P(0i + 2 ) = P(0i) + P(0 2 ). (36) 
Finally, for the case of N summands we obtain (compare with [TJ1 [151 fTT]): 

(jv \ iV 
^0„ + A' = ^£(0 n )+A, (37) 
ra=l / n=l 
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(N \ N 

n=l / n=l 

From formulas (I3T1) and (j3"8|) we see that the energy density and the pressure are sums 
of "individual" energy densities and pressures respectively and have no crossing term. 

In the case of an arbitrary metric g a p and a scalar field (f> n (t, xi, X2, £3), which satisfies 
eq. (JTT1) . we obtain that 



T^(0 n ) = -g a/3 I U 2 n - ^d^ n d^ n - Ue^^f ) - CU,, ,„ 



-c^ a/3 / rfp [(e( 1+ ") D ^ n )(n g e( 1 -") [:] ^ n ) + (^e( 1+ ") D ^ n )(^e^ D ^ n )] + 



1 

+ 2cJ dp{d a e^ u ^ n ){dpt 



The energy-momentum tensor, which corresponds to the function (jTUl) . is as follows 



N \ N 



T a/ 3 = T a(3 0n 1 + fl^A' = ^ T a/3 ((f) n ) + g a pA'. (39) 



K n=l / n=l 



4.2 Energy Density and Pressure for real a 

As we have seen in Section 3 for some values of parameters £ and c eq. (1121) has real roots. 
We denote as r\ m the value of r\ a for real a = m: 

2m? _ 

Vm = e(l + 2m 2 ) - 2 = (1 + 2m 2 ) - 2. (40) 

If and only if c > 1, then there exists the interval of < m 2 < vr? max ^ on which i] m < 
0. Some part of this interval is not physical, because g(m 2 , c) < on this part. The 
straightforward calculations (compare with [TJ]) show that at the point 

<a* = ~ \ ~ \W-X (-^) , (41) 

we obtain r\ m {ra max ) = 0. So, for c > 1 and £ 2 < £^ ax we have two positive roots of ( |T2|) : mi 
and m 2 > m\, with i] mi < and r] m2 > 0. In the next section we use this fact to construct 
a quintom local model with one tachyon real scalar field, which corresponds to r] m2 , and one 
phantom real scalar field, which corresponds to r] mi . For different values of c the function 
p m = m 2 r] m is presented in Fig. [2J 
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Figure 2: The dependence of p m on m at c = 1/2 (right), c = 1 (center) and c = 2 (left). 

5 Construction of solutions in the Friedmann— Robert- 
son— Walker metric 



5.1 Equations of motion and Friedmann equations 

In the spatially flat Friedmann-Robertson- Walker Universe we get the following equation of 
motion for the space homogeneous scalar field <fi 

{?V + l)e- 2V <P = c<p. (42) 



The Friedmann equations have the following form 



3H 2 = — E, 



mi 



H= - 



2m 2 



where dot denotes the time derivative (H = d H). 

The second equation of system (j4~3l) is the nonlinear integral equation in H(t): 



H= - 



m 



2 i 2 



(«9 O 0) 2 - c / (d e (1+ ^) (d Q e^cj>) dp 



(43) 



(44) 



Let us make an assumption, that <p(t) and H(t) satisfy the following equation 

V(j)= -a 2 (f), (45) 

where a is a root of eq. (TT21 . 

In this case eq. (jl2l is solved. Using formulas ff26l) and (T27j), we rewrite system (T43]) in 
the following form: 

(46) 



2m 2 , 



H= - 



2m 2 
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It is easy to check that (H5|) is a consequence of system (14"6T) . Instead of (HoT) we can consider 
the following third order system: 



+ 3#0 = a 2 0, 



# = - 



i2 



2m^ 



This system has the following integral of motion: 



h = 3H 2 



1m 2 v 



9 - a V 



(47) 



:A', 



mi 



(48 



therefore, choosing the initial date for ( 1471) one fixes the value of A'. 

So, our assumption allows to transform a system with a nonlocal scalar field into a system 
with a local one. In the same way we obtain systems with two or more local fields. Let 



N 



(49) 



n=l 



where all <j) n (t) are solutions of fT4"5l) with the same function H(t) and different values of a: 
a = a n . If all a n (n = 1..N) are different roots of (Tl2|) . then system (|4"3I) transforms into the 
following system with iV scalar fields: 



3H 2 



H 



i-ff 



2m » \ 

\n=l 

2m ? Its 



2 _ 2/2 



Va n <fil 



In the case of two real roots «i > and a%> ot\. 



+ 2A' 



(50) 



3H 2 



2mj 



1? - otic, 



+ 7], 



a ■> 



H= - 



2m* 



Va,4> 



i + Vc 



a 2/2 

\ ~ «202 



2A' 



(51) 



we have obtained that rj ai < and i] a2 > 0. Therefore the corresponding two-field model is a 
quintom one, in other words, includes one phantom scalar field (r] ai < ) and one scalar field 
with the canonical kinetic term (rj a2 > ) and with the tachyon mass term (a^c^ > 0)- The 
SFT inspired nonlinear local quintom models and their exact solutions have been studied, for 
example, in J37J[38]- To obtain exact solutions with physically important properties usually 
one should add some additional terms in the potential, which tend to zero in the limit of the 
flat space-time [151 E3 EEJ EH] . It is interesting that system (T46j) allows to find a physically 
important exact solution without adding any term in the potential. 
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5.2 Exact Solution in the case N = 1 



Let us consider system (j46J) with real a. Two exact nontrivial real solutions of this system 
have been presented in [16]. In our notations these solutions are the following: 



At a 7^ and r\ a < 



a 2 



<P(t) = A(t - t ), A' = - A , H(t) = —(t - t ), (52) 



where 



A = ±J- 2 ^, (53) 



3r L 



to is an arbitrary constant. 

At a = 0, A' = and r] a = £ 2 - 2 > 



'2m 2 I 
= ± v /_£ln(t_t ) +d, ff(t) = ^y, (54) 

where to an d are arbitrary constants. Note that the root a = exists if and only if 
c= 1. 

In this paper we present a new solution, which looks more realistic for the SFT inspired 
cosmological model. At present time one of the possible scenarios of the Universe evolution 
considers the Universe to be a D3-brane (3 spatial and one time variable) embedded in 
higher-dimensional space-time. This D-brane is unstable and does evolve to the stable state. 
This process is described by the dynamics of the open string, which ends are attached to 
the brane (see reviews [7] and references therein). A phantom scalar field is an open string 
theory tachyon. According to the Sen's conjecture [28] this tachyon describes brane decay, 
at which a slow transition in a stable vacuum takes place. This vacuum is characterized 
by the absence of open string states, i.e. corresponds to states of the closed string. This 
picture allows us to specify the asymptotic conditions for the scalar field. We assume that 
the phantom field <f>(t) smoothly rolls from the unstable perturbative vacuum (0 = 0) to a 
nonperturbative one, for example = A , where A is a nonzero constant, and stops there. 
It is easy to see that exact solutions, presented in [16j do not satisfy these conditions. 

At c = 1 our model is a nonlocal model for the dilaton coupling to the gravitation 
field. Its distinguished feature is the invariance under the shift of the dilaton field to a 
constant. In this case one of solutions of eq. ffl2|) is a = 0. Summing the first and the second 
equations of (JIB]) , we obtain: 

A' 

H = — - 3H 2 . (55) 
If A' > 0, then we obtain a real solution: 

#i(0 = l/^ tanh l J^(t-t ) ), (56) 
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where to is an arbitrary real constant. 

It is easy to see that Hi(t) > for any t, hence, from the second equation of fj46l) we 
obtain that <fi(t) can be real scalar field only if it is a phantom one (r] a < 0, that is equivalent 
to £ 2 < 2). The explicit form of <p(t) is as follows: 



<j> 1 {t) = ±\ 



■ arctan I sinh 



3(2 -e 



/3A' 

-5"(*-*o) 



+ Q 



mi 



2; 



(57) 



where is an arbitrary constant. Functions Hi(t) and 4>i(t) are presented in Fig. [3j 
1.0- 




-1.0-^ 



Figure 3: The functions Hi(i) (right) and 0i(t) (left) at A' 
C 2 = 0. 




-1.5— I 
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1, t, 



5 fc 



and 



(5f 



The Hubble parameter i?i(t) is a monotonically increasing function, so, using that 

2£T a 

we obtain u> < — 1. So, solution (1561) corresponds to phantom dark energy. Note that we 
have found two-parameter set of exact solutions at any A' > 0. In other words, at any 
A' > we have found the general solution of fj46|) . which correspond to a = 0. At A' = 
the solution fl56l) transforms to a constant. In the case A' = the general solution has been 
found in [T6] . 

In the case A' < we obtain the following general solution: 



H 2 (t) 



-A' 

3m^ 



tan 



3A' 



mi 



(t - 1 ) 



( 



± 



3(e - 2) 



- — -arctanh 



cos 



Co. 



(59) 



(60) 
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This solution is real at £ 2 > 2. It is interesting that the type of solutions essentially 
depends on sign of A'. The solution with the SFT inspired boundary conditions corresponds 
to A' > 0. 



6 Conclusions 

We have studied the SFT inspired linear nonlocal model. This model has an infinite number 
of higher derivative terms and are characterized by two positive parameters: £ 2 and c. 
For particular cases of the parameters £ 2 and c the corresponding actions describe linear 
approximations to either the bosonic or nonBPS fermionic cubic SFT as well as to the 
nonpolynomial SFT. 

Roots of the characteristic equation do not depend on the form of the metric and this 
property allows us to study properties of energy density and pressure. We have found that in 
an arbitrary metric the energy-momentum tensor for an arbitrary N-mode solution is a sum of 
the energy-momentum tensors for the corresponding one-mode solutions. In the Friedmann- 
Robert son- Walker spatially flat metric the pressure for a one-mode solution corresponding 
to a real root can be positive or negative, depending on parameters of our nonlocal model. 
Namely, for c < 1 the one mode pressure is positive and for c > 1 it could be negative or 
positive. 

The investigation performed in this paper shows that the general field equations in linear 
nonlocal models admit an equivalent description in terms of local theory and as a conse- 
quences we have representations (37) and (38) for the energy and pressure. This calculation 
also supports the use of the Ostrogradski representation for our system in the case of arbi- 
trary metric. 

In distinguish to our previous paper [J5] we do not use any approximation scheme and do 
not add any terms in the potential. We have shown that our linear model with one nonlocal 
scalar field generates an infinite number of local models. These models can be studied 
numerically and we plan to present this analysis in future papers. Some of these models 
have been solved explicitly and, hence, special exact solutions for nonlocal model in the 
Friedmann-Robertson- Walker metric have been obtained. In particular we have constructed 
an exact kink-like solution, which correspond to monotonically increasing Universe with 
phantom dark energy. Note that the obtained behaviour of the Hubble parameter is close 
to behavior of the Hubble parameter in the nonlinear nonlocal model [8] , which recently has 
been obtained numerically [19|. 
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